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Abstract 
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Extracellular local field potentials (LFP) are usually modeled as arising from a set of current sources 
embedded in a homogeneous extracellular medium. Although this formalism can successfully model 
several properties of LFPs, it does not account for their frequency-dependent attenuation with dis- 
tance, a property essential to correctly model extracellular spikes. Here we derive expressions for 
the extracellular potential that include this frequency-dependent attenuation. We first show that, if 
the extracellular conductivity is non-homogeneous, there is induction of non-homogeneous charge 
densities which may result in a low-pass filter. We next derive a simplified model consisting of a 
punctual (or spherical) current source with spherically- symmetric conductivity /permittivity gradients 
around the source. We analyze the effect of different radial profiles of conductivity and permittivity 
on the frequency-filtering behavior of this model. We show that this simple model generally displays 
low-pass filtering behavior, in which fast electrical events (such as Na + -mediated action potentials) 
attenuate very steeply with distance, while slower (K + -mediated) events propagate over larger dis- 
tances in extracellular space, in qualitative agreement with experimental observations. This simple 
model can be used to obtain frequency-dependent extracellular field potentials without taking into 
account explicitly the complex folding of extracellular space. 
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Extracellular potentials, such as local field potentials (LFPs) or the electroencephalogram (EEG), are 
routinely measured in electrophysiological experiments. The fact that action potentials have a limited 
participation to the genesis of the EEG or LFPs was noted from early studies. Bremer (1938, 1949) 
proposed that the EEG is generated by non-propagating potentials, based on the mismatch of time 
course between EEG waves and action potentials. Eccles (195 1) proposed that LFP and EEG activities 
are generated by summated postsynaptic potentials arising from the synchronized excitation of cortical 
neurons. Intracellular recordings from cortical neurons later demonstrated a close correspondence 
between EEG/LFP activity and synaptic potentials (Klee et al., 1965; Creutzfeldt et al., 1996a, 1996b). 
The current view is that EEG and LFPs are generated by synchronized synaptic currents arising on 
cortical neurons, possibly through the formation of dipoles (Nunez, 1981; Niedermeyer and Lopes da 
Silva, 1998). 

The fact that action potentials do not participate to EEG-related activities indicate strong frequency- 
filtering properties of cortical tissue. High frequencies (greater than « 100 Hz), such as that produced 
by action potentials, are subject to a severe attenuation, and therefore are visible only for electrodes 
immediately adjacent to the recorded cell. On the other hand, low-frequency events, such as synaptic 
potentials, attenuate less with distance. These events can therefore propagate over large distances in 
extracellular space and be recordable as far as on the surface of the scalp, where they can participate 
in the genesis of the EEG. This frequency-dependent behavior is also seen routinely in extracellular 
unit recordings: the amplitude of extracellularly-recorded spikes is very sensitive to the position of 
the electrode, but slow events show much less sensitivity to the position. In other words, an extracel- 
lular electrode records slow events that originate from a large number of neighboring neurons, while 
the action potentials are recorded only for the cell(s) immediately adjacent to the electrode. This 
fundamental property allows to resolve single units from extracellular recordings. 

However, little is known about the physical basis of the frequency-dependent attenuation of extra- 
cellular potentials in cortex. By contrast to intracellular events, which biophysical mechanisms have 
been remarkably well characterized during the last 50 years (reviewed in Koch, 1999), comparatively 
little has been done to investigate the biophysical mechanisms underlying the genesis of extracellular 
field potentials (see review by Nunez, 1981). The reason is that LFPs result from complex interac- 
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tions involving many factors, such as the spatial distribution of current sources, the spatial distribution 
of positive and negative electric charges (forming dipoles), their time evolution (dynamics), as well 
as the conductive and permittivity properties of the extracellular medium. One of the simplest and 
widely used model of LFP activity considers current sources embedded in a homogeneous extracel- 
lular medium (Nunez, 1981; Koch and Segev, 1998). Although this formalism has been successful in 
many instances (Rail and Shepherd, 1968; Klee and Rail, 1977; Protopapas et al., 1998; Destexhe, 
1998), it does not account for the frequency-dependent attenuation and therefore is inadequate for 
modeling extracellular field potentials including spike activity. 

In this paper, we would like to investigate possible physical grounds for the frequency-filtering 
properties of LFPs. We start from first principles (Maxwell equations) and consider different condi- 
tions of current sources and extracellular media. We delineate the cases leading to frequency-filtering 
properties consistent with physiological data. We show that the assumption of a resistive homoge- 
neous extracellular medium cannot account for the frequency-dependent attenuation. It is necessary 
to take into account the inhomogeneous structure of the extracellular medium (in both permittivity and 
conductivity) in order to account for frequency-dependent attenuation. We next analyze a simplified 
representation of current sources in non-homogeneous media, and provide a simplified model which 
could be applied to simulate extracellular field potentials without using complex representations of 
extracellular space. We terminate by showing a concrete example of the genesis of extracellular LFPs 
from a conductance-based spiking neuron model. 

MATERIAL AND METHODS 

We will first develop a general formalism to express the temporal variations of extracellular potential, 
as well as a simple model in which most of the calculations can be done analytically (see Section 
General theory; see also Appendix 1 and 2 for details). We will next explore this simplified model 
numerically to illustrate its frequency-filtering behavior (Section Numerical simulations), in which 
we have performed two computations: (i) Calculate the impedance: the impedance is given by an 
integral, which was evaluated numerically by standard numerical integration routines, (ii) Convert 
time-dependent functions into frequency-spectra. These conversions were done via Fourier transfor- 
mation (as well as its reverse transformation), which were carried out in C using standard numerical 
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routines (Press et al., 1986). 

To test this formalism, we also considered a simple biophysical model of a spiking neuron con- 
taining voltage-dependent and synaptic conductances (last part of Section Numerical simulations). 
A single-compartment neuron was constructed and included conductance-based models of voltage- 
dependent conductances and synaptic conductances. This model was described by the following 
membrane equation: 

dV 

C '"^7 = -8l(V-E l ) -g Na (y-E Na ) -g Kd (V-E K ) - gM (V-E K ) - ge (V-E e ) , (1) 
dt 

where C m = 1 /jF/cm 2 is the specific membrane capacitance, gi = 4.52 x 10~ 5 S/cm 2 and Ei = - 
70 mV are the leak conductance and reversal potential. g^ a = 0.05 S/cm 2 and g^d = 0.01 S/cm 2 
are the voltage-dependent Na + and K + conductances responsible for action potentials and were de- 
scribed by a modified version of the Hodgkin & Huxley (1952) model. gM = 5 x 10~ 4 S/cm 2 is a 
slow voltage-dependent K + conductance responsible for spike-frequency adaptation. g e = 0.4 /jS is 
a fast glutamatergic (excitatory) synaptic conductance. The voltage-dependent conductances were 
described by conventional Hodgkin-Huxley type models adapted for modeling neocortical neurons, 
and the synaptic conductance was described by a first-order kinetic model of neurotransmitter binding 
to postsynaptic receptors. These models and their kinetic parameters were described in detail in a 
previous publication (Destexhe and Pare, 1999). All numerical simulations were performed using the 
NEURON simulation environment (Hines and Carnevale, 1997). 



GENERAL THEORY 

In this section, we outline the main features of the model starting from first principles (Maxwell 
equations). We will consider a number of different special cases and derive a simplified model with 
radial (spherical) symmetry. In Section Numerical simulations, we will investigate numerically the 
behavior of this simplified model. 

In Maxwell's theory, the electric properties of a conductive medium are determined by two pa- 
rameters, conductivity a and permittivity £. While conductivity quantifies the local relation between 
the electric field and the current, permittivity characterizes the response of the system in terms of 
separation of opposite charges (polarisation) in the presence of an electric field. Maxwell's theory of 
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electromagnetism allows one to compute electric and magnetic fields or potentials, for a given distri- 
bution of charges and electric currents. Because charges move very slowly in biological media, the 
effects of magnetic fields are very small compared to that of the electric field and will be neglected 
here. 

One of Maxwell's equations is Gauss' law. 

V-(£E)=p. (2) 

Here E denotes is the electric field, D = EE is the called displacement and p is the charge density of 
the extracellular medium, also allowed to vary slowly in time. 

The continuity equation relates the current density j to the charge density p: 

V.J + |=0 (3) 

which states a balance between the electric flux into some volume and the change of the total charge 
in this volume. In other words, no charge will get lost. 
Finally, there is Ohm's law: 

j = oE, (4) 

which states that the relation between the local electric field and the local current density can be 
described by a single macroscopic parameter a. Here, we assume that the electric field and the current 
are parallel locally (i.e., at any given point in the medium). The coefficient relating E and j, the 
conductivity a, is a scalar, which is justified by macroscopic measurements (Ranck, 1963). 
Combining Eqs.|3] and |4] yields 

V-(aE) + ^ = 0. (5) 

This equation describes charge conservation in differential form. 

In the following, we assume that electric currents are distributed on the surface of the membrane 
(ionic currents), and that these currents are allowed to vary in time. We introduce the the electric 
potential V (also called extracellular potential) which obeys 

E = - W . (6) 

Below, we successively consider different cases of increasing complexity, starting with a homoge- 
neous extracellular medium, then going over to non-homogeneous media. 
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Homogeneous extracellular medium 

Consider a membrane embedded in a homogeneous extracellular medium, with conductivity o and 
permittivity £ being held constant in space and time. As shown in Appendix 3, using the assumption 
a = Go = const., £ = £o = const., we get for each spectral component 

+ AV. = AV.=0. (7) 



Hence we find 



AV«, = -^ = 0. (8) 



Using the inverse Fourier transform, this yields 

AV(x,*) = -~ = 0. (9) 

One observes that the charge density p vanishes at the exterior of the sources. The solution depends 
on the geometry considered, its symmetries and boundary conditions. We consider different cases 
below. 

Spherical membranes 

As a particular case, let us consider a spherical membrane obeying the conditions outlined above 
(homogeneous medium of constant conductivity and permittivity), in which the potential is allowed 
to vary as a function of time. There are two classic cases of boundary conditions, which yield unique 
solutions of the Laplace equation: (i) Specification of the potential on a surface (Dirichlet boundary 
condition) and (ii) Specification of the derivative of the potential normal to the surface (Neumann 
boundary condition). The latter is equivalent to specifying the normal component of the electric field. 
Due to Ohm's law this means to specify the current density on the surface. In neurons, we consider 
the Neumann conditions in which current sources represent the ionic currents in the membrane. In 
this case, we have the following boundary conditions: 

j(r = R) = j n = oE(r = R) , (10) 

Here, we have assumed that the membrane is spherical of radius R (which is also equivalent to consider 
a point current source). The current density on the surface of the sphere is is proportional to the unit 
vector n normal to the surface of the sphere. 
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To analyze the behavior of the extracellular potential as a function of frequency, we perform a 
Fourier transformation with respect to time of the potential, 

/CO 
dtV(x,t)e iat . (11) 
-oo 

The potential V satisfies the Laplace equation (11). Because the Fourier transform is linear, each 
Fourier component V m also satisfies (10). The solution, satisfying boundary conditions (Eq. ITOt. is 
given by 

"•M = 5b- (12) 

where I® = ju>o/(4KR 2 ) is the total current for each frequency component. This shows that the po- 
tential is the same for all frequency components. Therefore, there is no frequency-dependence in this 
case. 

Thus, the general expression for field potentials resulting from a set of N current sources {/,•} of 
spherical symmetry is: 

1 N I 
47UG £ I*"*/ 

where x is a point in extracellular space, Xj is the location of the jth current source, and |x — x 7 | is 
the distance between x and Xj. This expression is widely used to model extracellular field potentials 
(Nunez, 1981; Koch and Segev, 1998). 

Cylindric membranes 

Because neuronal processes (dendrites, axons) are closer to cylinders, we considered cylindric mem- 
branes as a second particular case. Here, the procedure is similar to the above, but the geometry and 
symmetries are different. If we assume that the membrane potential respects cylindric symmetry (ie, 
V does not depend on the rotation angle around the cylinder axis), and is uniform on the surface of the 
membrane (isopotential compartment), then Eq.[9]can be solved for each component of the frequency 
spectrum and has a unique solution: 

V M (r) = V a (R) + j m (R) | Hr/R) , (14) 

where r is the distance perpendicular to the cylinder axis, while j a (R) and V a (R) are the (O-frequency 
components of the current density and potential at the surface of the cylinder. 
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In this case again, the extracellular potential is independent of the frequency. The same conclusion 
applies to membranes of arbitrary geometries taken in the same conditions because the (O-frequency 
component of the current density and potential obey Laplace equation. There is therefore no frequency 
dependence arising from single current sources in homogeneous media. 

Non-homogeneous extracellular medium 

We have shown above that, in a homogeneous medium with a current source of spherical or cylindric 
symmetry, the extracellular electric potential is the same for all frequency components, and there- 
fore cannot display frequency-dependent properties. We now turn to a possible source of frequency- 
dependent attenuation, namely the presence of inhomogeneities in the conductivity of the extracellular 
medium. 

Stationary currents in spherically-symmetric non-homogeneous medium 

Before investigating the general case, let us first consider the case of a static spherical current source 
embedded in a medium where the conductivity a conserves spherical symmetry, but varies as a func- 
tion of distance r (as above we assume that a does not depend on time). We also continue to assume 
that permittivity £ is homogeneous. If the total current flowing through the sphere of radius R is 
denoted by /, then the radial dependence of the current density is given by 

In this case, the charge density p is non-zero and is given by 

p = --j-Vlogo. (16) 
a 

Then Ohm's law implies for the spherically symmetric electric field 



/ 

4nr 2 o(r) 



E M = 1 2 — ~ *r ■ (17) 



The spherically symmetric electric potential is the obtained by integrating the electric field, giving 



V(r) = -l'dr'E(r')= /V ' . (18) 
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Details of the calculation can be found in Appendix 1 . 

This equation shows that the potential may decrease or even increase, depending on the spatial 
variations of a. An important consequence is that such net charge creates its own electric field (so- 
called secondary field), which will be analyzed in more detail below. What is the physical origin of 
this non-zero net charge? The current density behaves like J °= l/r 2 , (Eq. fT5t . and the electric field 
like £« 1/ (a(r)r 2 ), (Eq.fTTl. Consequently, there will be accumulations of charges in some regions 
of lower conductivity, similar to traffic jams. Consider a more realistic case in which the conductivity 
of the extracellular space is constant on average, but displays spatial fluctuations around this average. 
This could correspond for example to different processes and obstacles in the extracellular medium. 
In this case, the electric field , going like 1 /r 2 on average, fluctuates locally. This creates local areas 
of positive and negative charge, i.e. electric dipoles. Those dipoles also create a secondary electric 
field. However, to account for frequency dependence, time-varying current sources must necessarily 
be considered, in which case the situation is more complex, as analyzed in the next section. 

Time-varying currents in non-homogeneous medium 

Let us now consider the general case where both e and a are non-homogeneous in space, but constant 
in time. We assume that the current source is allowed to vary in time. The continuity equation implies 
that the charge density is also time-dependent. Ohm's law implies that the electric field has a time- 
dependence as well, and so will also the extracellular potential. Due to the inhomogeneity of £, the 
extracellular potential does no longer satisfy Poisson's equation. To study the frequency-dependence 
of the extracellular potential, we perform a Fourier transform of the electric field, the potential and 
likewise of the charge density p. 

pco, the component of frequency (0 of the temporal Fourier transform of the charge density p, 
satisfies 



This equation expresses the differential law of charge conservation for a given Fourier component. 
Now we consider the Gauss' law (Eq.EJ), the law of charge conservation in differential form (Eq. O, 
and carry out the Fourier transform with respect to time. Taking into account Eq.fT9lvields 




imp a . 



(19) 



(Wco)-(V(G + Koe)) 



(W ffl )-(Vlog(o + i(oe)). 



(20) 



O + ?G)£ 
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For details, see Appendix 2. 

This equation is general and applies to any particular symmetries (under the assumption of scalar 
conductivity). We consider below a series of special cases, as well as special symmetries. 

Special cases 

As a first special case, consider Eql20lwhen permittivity is constant. Then Ohm's law (Eq.|4|) implies 

AVco = -— • (21) 
£ 

Constant permittivity also implies V(o + zd)£) = V(a). Hence Eq.|20]takes the form 

A y. = _ W -- V ° = -E-. (22) 
O + Z(D£ £ 

We therefore observe the occurrence of a (complex) phase difference between the induced charge 
density p/£ and the current density j (recall: — VVco = Eco = jco/o). This effect depends on the 
frequency 0) of the Fourier component. Such phenomenon is well known from electric circuits of the 
RC type, where in general a phase difference between potential and current is observed. In particular, 
if the potential vanishes at some time t, the electric charge density will not immediately go to zero. 

Eq.[22]shows that for high enough frequency the induced charge density goes to zero. On the other 
hand, for low-frequency phenomena, the charge density will carry out large fluctuations and will be 
sensitive to spatial fluctuations of conductivity. This has important consequences for interpreting LFP 
activity (see Section Frequency-filtering properties of non-homogeneous media). 

As another special case, consider Eq. |20] when the conductivity is constant. The law of charge 
conservation in differential form Eq.|5]then becomes 

AVco = -^- (23) 
a dt 

Constant conductivity also implies V(o + z(Q£) = V(z'o)£). Hence Eq.|20|takes the form 

_vvvv^ = iap £ = _ iv 

O + Z0)£ O dt O 
This means in the limit of low frequencies that there are no current sinks or sources. 

A third noteworthy special case is when both permittivity and conductivity are non-homogeneous, 
but have a fixed ratio: 

— = const. (25) 
a 
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Under those circumstances one obtains 



(VVco) • (Vc) 



(W ra ) ■ (Vlog(a)) . 



(26) 



co — 



c 



All frequency-dependence cancels out, the potential becomes frequency independent. 

Frequency-filtering properties of non-homogeneous media 

We now analyze the frequency-filtering properties of Eq. |20| Consider this equation in two limit 
cases: (i) eo<< c/e (low frequency limit); in this case, log (a + ia>e) « log (a), the solution becomes 
independent of the permittivity e and is determined only by the conductivity c. (ii) CO >> c/e (high 
frequency limit); in this case, log(c-Mcoe) ~ log(zeoe), and permittivity only determines the solution. 
These two cases will be considered in more detail below. The critical frequency around which this 
transition will occur depends on the relative values of c and e, and one can define the following critical 
frequency f cr 



As an example, consider the value of average resistivity p res (inverse of conductivity c) measured 
in rabbit cerebral cortex (Ranck, 1963), giving p res = 3£l m. Taking the permittivity of salt water 
(e = 7 x Kr 10 F/m), gives a critical frequency of f cr of about 10 10 Hz. Thus, this analysis shows that 
the behavior will be similar for low and high frequency limits. However, if one evaluates 2no/e for 
a resting membrane (closed ion channels; p ~ l0 9 Q.m and e ~ 10~ 10 F/m), one finds for the critical 
frequency a value in the range between and 100 Hz. The phenomenon of induced charges will be 
likely to play a role in the frequency range of synaptic inputs in cerebral cortex (0 to 40 Hz). On the 
other hand, higher frequencies (> 100 Hz) - such as action potentials - are likely to cause negligible 
variations in charge density. 

Further, we can compare Eq. [20] to cases where there is no frequency dependence (i.e., Laplace 
Eq. At the limit of high frequencies, the left term of Eq. ESI vanishes and this equation becomes 
equivalent to Eq.|5J showing that for high frequencies, one recovers the same behavior as for a homo- 
geneous medium. To have a low-pass filter, similar to what is observed from extracellular recordings, 
one must have a situation in which the attenuation of the potential at low frequencies must be less than 
for homogeneous media. Inspection of the left term of Eq.|5D]shows that the attenuation can be either 



27100, 



a 



(27) 
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less or more pronounced, resulting in low- or high-pass filters. The type of filter will depend on the 
behavior of the gradients of conductivity and permittivity. This behavior will be analyzed numerically 
in more detail later (see Section Numerical simulations). 

Time-varying currents in spherically-symmetric non-homogeneous medium 

In order to calculate the extracellular potential V generated by time-varying currents in non-homo- 
geneous media, Eq. must be integrated by incorporating details about the particular geometry 
of current sources and extracellular properties. A general method for solving this problem is, e.g., 
finite-element analysis, which allows to explicitly incorporate the complex shape and composition 
of extracellular space around neurons. However, this approach requires to integrate complex mor- 
phological data and appropriate simulation tools. We defer this to a future study. In order to have a 
model of LFPs applicable to standard neuron models, we follow here a simpler approach, based on 
the following simplification: we consider that the variations of conductivity and permittivity have a 
radial symmetry in the vicinity of the current sources. This simplification allows us to obtain sim- 
pler expressions of the extracellular potential, still displaying frequency dependence, and apply this 
formalism using standard simulation tools. 

Consider Eq.ESlfor the case of a spherically-symmetric system. Then the potential obeys 

d 2 V a 2dV m 1 d(o + /coe) dV a ^ 

, 9 H : \~-, — r : — — U. (ZQJ 

dr L r dr (o + i(D£) dr dr 
Integrating this equation gives the following relation between two points r\ and r2 in the extracellular 

space, 

A ^(l) [o(n) + ime(n)] = r\ ^(r 2 ) [o(r 2 ) + toe(r 2 )] . (29) 
dr dr 



This can be verified by differentiating this equation with respect to r, which then yields Eq. |28J Inte- 
grating Eq.E3once more yields 

v-(n)-v^ + g W £ ^;:a - 

This can be seen most easily by differentiating Eq. |30|to yield Eq. |29j In particular, if r 2 = R, — ^ (r 2 ) 
represents the electric field at the surface of the sphere of radius R, which by Ohm's law is related to 
the current density at R, we obtain 
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If we assume that the extracellular potential vanishes at large distances (Vco(°°) = 0), we have: 

V(o{°°) = = Vco(a — — - / dr — T — - — — — . (32) 

K ' K ' 4no(R)J R r' 2 cj(r') + i(D e(r') 

which allows to eliminate V a (R) from Eq.|3H leading to: 

This will be the main equation that forms the basis of our simplified model of LFP. Solving this equa- 
tion for o and £ constant leads to the expression found above (Eq. IT"2"1) for homogeneous extracellular 
media. In the numerical part (Section Numerical simulations), we will solve this equation for different 
spatial profiles of o and 8. To this end, it is useful to define the impedance: 

~ / n 1 r , / 1 o(R) + z'G) e(R) 

Z » r i = ~j T^r / dr' \ )-( . 34 

Then Eq. ^becomes 

V a (n) = Z ra (ri) 7 ro . (35) 

The impedance is therefore the "filter" applied to the co-frequency component of the current source, 
to yield the corresponding frequency component of the extracellular potential. In the next section, we 
will examine the frequency-filtering properties of different extracellular media by calculating numer- 
ically the impedance for different cases of spatial inhomogeneities of conductivity and permittivity. 



NUMERICAL SIMULATIONS 

In this section, we use the expressions of the extracellular potential obtained above. In particular 
we analyze the behavior of the extracellular potential generated by a current source in a spherically- 
symmetric non-homogeneous medium (Eq. l33l). and its associated impedance (Eq. l34l). We investi- 
gated the frequency-filtering properties obtained for different cases of increasing complexity of the 
radial profile of a and £. The goal is to determine the conditions of spatial variations of conductiv- 
ity and permittivity for which the frequency-filtering properties are consistent with physiological data. 
We terminate by an application of this model to calculating the LFP generated by a conductance-based 
spiking neuron model. 
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Parameters 

Precise experimental data on the variations of permittivity £ and conductivity o in the extracellular 
medium have not been measured so far. However, averaged values of these parameters are available 
from macroscopic measurements. A value for a, averaged over large extracellular distances, o av , 
was measured by Ranck (1963) and was between 0.28 S/m and 0.43 S/m, for 5 Hz and 5 kHz, 
respectively. The macroscopic frequency dependence of conductivity seems therefore relatively weak. 
However, the situation is different microscopically. As reviewed in Nunez (1981), the conductivity of 
the CSF fluid is 1.56 S/m while the typical conductivity of membranes is 3.5 x 10~ 9 S/m. This value 
was obtained from the resting (leak) membrane conductance of cortical neurons, typically around 
4.5 x 10 5 S/cm 2 , multiplied by the thickness of the membrane (7-8 nm; Peters et al., 1991). At 
microscopic scales, there is therefore approximately 9 orders of magnitude variations of conductivity. 

Permittivity variations are not so dramatic. Fluids have higher permittivity, for example it is 
about 7 x 10~ 10 F jm for sea water. Membranes have a permittivity of about 7.5 x 10 -11 F/m. The 
latter value was derived from the specific capacitance of membranes, C = 1 /j F/cm 2 (Johnston and 
Wu, 1997), and assuming a membrane thickness of 7.5 nm (Peters et al., 1991). Because those 
variations are small compared to the variations of conductivity, it is a good approximation to consider 
the permittivity as a constant. In the following, we will use the reference value of £ = 10~ 10 F/m. 

In the following we will use normalized values for conductivity c(r) /c(R) and permittivity e(r) /c(R) . 
Because the membrane is always surrounded by extracellular fluid, o(R) = 1.56 S/m and the normal- 
ized conductivity o(r)/o(R) therefore varies between 1 and about 2 x 10~ 9 . Similarly, the normalized 
(constant) value of permittivity will be e(r)/a(R) = 6 x 10~ u s. 

Frequency-filtering properties of spherically-symmetric media 

We calculated numerically the impedance (Eq.l34l) for different cases of spatial variations of conduc- 
tivity and permittivity. In all cases we assumed a current source with spherical geometry, characterized 
by radius R, and that o and £ vary according to a radial (spherical) symmetry around this current source 
(see scheme in Fig. UK). For each case, we represented the normalized impedance 



2 a {r)=Z at {r)/Z at {R) , 



(36) 
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which allows better comparison between the different cases. Because the value of impedance does not 
depend on the absolute value of permittivity and conductivity (dividing o and £ by a constant factor 
does not change Z m (r) in Eq.O?t. we used the normalized conductivity a(r)/a(7?) and the normalized 
permittivity e(r)/a(2?) defined above. It is also convenient to represent all distances in units of R, 
although we also considered absolute values of distances (see below). 

We first investigated a simple case of smooth variations of those parameters, to illustrate the dif- 
ferent types of frequency filtering that can be obtained in this model. The profiles of conductivity 
and permittivity are shown in Figs.[Tj3 and C. These curves tend to the same asymptotic value for 
large distances. The corresponding impedance is shown as a function of frequency / in Figs.[TJ)-F 
(see Appendix 3 for details of the method). When the ratio o/£ is kept constant (Fig.[ll)-F, dotted 
line), there is no frequency dependence as analyzed above in Section Special cases. In the case of 
a decreasing conductivity with distance combined with constant permittivity, one has a high-pass fil- 
ter (Fig. UJD-F, dashed line). By contrast, a low-pass filter is observed if an increasing conductivity 
with distance is combined with a constant permittivity (Fig. UJD-F, solid line). Thus, there is a clear 
frequency-dependent behavior when a and/or £ vary as function of distance r, if the ratio o/e is not 
constant. This also shows that low- and high-pass filters are both possible, depending on the exact 
form of the function o(r) and e(r). The impedance can also have a non-zero imaginary part, which 
means that beyond resistivity, the medium has is also capacitive properties. In this case, there will be 
a phase difference between the potential and the current. 

We next considered a case characterized by a localized drop of conductivity (Fig.EK) while per- 
mittivity was kept constant (Fig. E^). The resulting impedance measured at different distances from 
the source is shown in Fig.|2t-E as a function of frequency /. In this case, for distances around the 
conductivity drop, there is a moderate frequency dependence with low-pass characteristics (Fig.Et- 
E, dotted and dashed lines). However, for larger distances, the imaginary part is zero and there is no 
frequency dependence (Fig.|2l!-E, solid lines). This is explained by the fact that for large distances 
a(r) = o(R) and e(r) = e(R). This behavior can also be seen in the attenuation of the different fre- 
quency components illustrated in Fig. . There is a different attenuation only for distances around 
the region where conductivity varies. 

Because the extracellular space is composed of alternating fluids and membranes (Peters et al., 
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1991), which have high and low conductivity, respectively, we have next considered the situation 
where conductivity fluctuates periodically with distance (Fig. [3l> . Considering a cosine function 
of conductivity (Fig. QK) with constant permittivity (Fig. |3f3) leads to a rather strong frequency- 
dependent attenuation (Fig.|3l!-E) with low-pass characteristics. There was a strong attenuation with 
distance for all frequencies (Fig. Of). Very similar results were obtained with other periodic functions 
(for example by replacing cos by sin in the function used in Fig.|3J\), different oscillation periods, or 
even for damped oscillations of conductivity (not shown). 

It could be argued that although fluids and membranes alternate in extracellular space, there is an 
efficient diffusion of ions only in the extracellular fluid around the membrane. For larger distances, 
diffusion becomes increasingly difficult because of the increased probability of meeting obstacles. In 
this case, conductivity would be highest around the source and progressively decrease to an "average" 
conductivity level for larger distances. This situation is illustrated in Fig.HJ We have considered that 
the conductivity is highest at the source, then decreases exponentially with distance with a space con- 
stant X (Fig.|4|\; note that in this case, real distances were used). Permittivity was constant (Fig.|4(3). 
The resulting impedance displayed pronounced frequency-filtering properties with low-pass charac- 
teristics (Fig.Hp-E). In particular, the attenuation with distance revealed strong differences between 
low and high frequencies of the spectrum (Fig.|4f). Similar results can be obtained with other de- 
creasing functions of connectivity (not shown). 

The above examples show that there can be a strong frequency-filtering behavior, with low-pass 
characteristics as observed in experiments. However, although these examples show a more effective 
filtering for high frequencies, it still remains to be shown that the high frequencies attenuate more 
steeply with distance compared to low frequencies. To this end, we define the quantity: 

Qioo = Z m (r)/Zi(r) , (37) 

where Z\ and Z\qq are the impedances computed at 1 Hz and 100 Hz, respectively. This ratio quanti- 
fies the differential filtering of fast and slow frequencies as a function of distance r. Fig.|5j\ displays 
the <2ioo values obtained for some of the examples considered above. In the case of a localized drop 
of conductivity (Fig.|3\, Drop), there was an effect of distance for r < \6R, then the <2ioo remained 
equal to unity for further distances. This behavior is in agreement with the impedance shown in Fig. El 
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in which case there was no frequency filtering for r > \6R. For oscillatory conductivities (Fig. [5]\, 
Osc), the <2ioo was always < 1, consistent with the low-pass frequency-filtering behavior observed in 
Fig.0 However, the <2ioo oscillated around a value of 0.6 and did not further decrease with distance. 
Thus, in this case, although there was a clear low-pass filtering behavior, all frequencies still contribute 
by the same relative amount to the extracellular potential, regardless of distance. On the other hand, 
with exponential decay of conductivity, the <2ioo monotonically decreased with distance (Fig. |5J\, 
Exp). Thus, this case shows both low-pass filtering behavior (Fig. |U) and a stronger attenuation of 
high frequencies compared to low frequencies (Fig. [5j\, Exp), which is in qualitative agreement with 
experiments. Analyzing exponentially-decaying conductivities of different space constants (Fig.|5j3) 
revealed that the various patterns of distance dependence approximately followed the pattern of con- 
ductivity (Fig.|5£). This type of conductivity profile is relatively simple and plausible, and will be the 
one considered in the biophysical model investigated below. 

Biophysical model of the frequency-filtering properties of local field potentials 

We have applied the above formalism to model the frequency dependence of the extracellular field 
potentials stemming from a conductance-based spiking neuron model. The details about the model are 
given in the Material and methods section, while the details of the calculation of the extracellular LFP 
(general for any current source) is given in Appendix 3. The profile of conductivity and permittivity 
used is that of Fig. |4j We calculated the total membrane current generated by a single-compartment 
model of an adapting cortical neuron, containing voltage-dependent Na + and K + conductances for 
generating action potentials and a slow voltage-dependent K + conductance responsible for spike- 
frequency adaptation. The model also contained a fast glutamatergic excitatory synaptic conductance, 
which was adjusted to evoke a post-synaptic potential just above threshold, in order to evoke a single 
action potential (Fig. 0A)- The total membrane current (Fig. m$) was calculated and stored in order 
to calculate its Fourier transform (power spectral density shown in Fig. 03). The impedance of the 
extracellular medium (Fig. (6jD) was calculated using absolute values of the parameters (Eq.l34l). 

This model was used to calculate the field potentials at different radial distances assuming the 
neuron was a spherical source (radius of 105 /mi). The extracellular potential is indicated for 5, 100, 
500 and 1000 fim away from the source (see Fig. [6^) and strong frequency filtering properties are 
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apparent: the fast negative deflection of extracellular voltage showed a steep attenuation and almost 
disappeared at 1000 /urn (although it had the highest amplitude at 5 /urn). In contrast, the slow positive 
deflection of the extracellular potential showed less attenuation with distance and became dominant 
at large distances (500 and 1000 in the example of Fig.|(Ji). 

Thus, this simple example illustrates that the approach provided here can lead to a relatively simple 
model to calculate local field potentials with frequency filtering properties. The exact profiles of 
filtering and attenuation depend on the exact shape of the gradients of conductivity/permittivity as 
well as on the spherical symmetry inherent to this model. 

DISCUSSION 

In this paper, we have provided a model of extracellular field potentials in non-homogeneous media. 
We discuss here the validity of this model, how it relates to previous studies, and what perspectives 
are provided. 

The theoretical analysis outlined in Section General theory shows that inhomogeneities of ex- 
tracellular space (with respect to conductivity and/or permittivity) is a possible cause for frequency 
filtering. In general, non-homogeneous extracellular media will differently affect the attenuation of 
the various frequency components of the current sources, and can lead to high-pass or low-pass filters 
depending on the gradients of conductivity and permittivity. The composition of extracellular space 
is made from the alternance of fluids and membranes (Peters et al., 1991). Because these media have 
very different conductivity and permittivity, one may expect that the extracellular space is necessarily 
highly non-homogeneous. Therefore, the structural composition of extracellular space is very likely 
to be a main determinant of the frequency-filtering properties of LFPs. In addition, the conductivity of 
the extracellular fluid directly beneath the membrane depends on the ionic concentrations present. It 
turns out that the extracellular ionic concentrations may vary in time, in an activity-dependent manner 
(reviewed in Amzica, 2002). Therefore, it is also likely that there is an activity dependent contribution 
to the filtering properties of the extracellular medium. Here, we did not consider such time-dependent 
variations of conductivity, but this type of contribution is certainly worth to be considered by future 
theoretical work. 

To correctly simulate the frequency-filtering behavior due to extracellular inhomogeneity, the ex- 
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tracellular potential should be calculated by a model incorporating details about the three-dimensional 
composition of the extracellular medium. Such type of simulations should use methods such as finite- 
element analysis. However, the complexity of this type of analysis, and of the data it requires, makes 
such simulations inaccessible to standard models. In addition, this requires orders of magnitude differ- 
ences in computational power needs. For these reasons, we have considered the option of generating 
a simplified model under some approximation. We assumed that the geometry of extracellular inho- 
mogeneities is spheric around the current source. In this condition, one can obtain relatively simple 
expressions of the extracellular potential such as Eq. |33] Not only this expression is amenable to 
theoretical analysis, but it is also sufficiently simple to be applied to current neuron models which do 
not have an explicit representation of extracellular space. 

The drawback of this method is that it considers an un-realistic (radial) distribution of inhomo- 
geneities in extracellular space, which will necessarily affect the frequency-filtering properties pro- 
duced by the model. However, it should be possible to calculate the "average" radial variations of 
conductivity and permittivity by averaging the profiles of a and £ in all directions emanating from neu- 
ronal membranes using three-dimensional reconstructions of the neuropil. Another direction would 
be to measure experimentally the profiles and distributions of a and £, but such data are not cur- 
rently available. Here, we have used "heuristic" profiles of conductivity and permittivity, which gives 
frequency-filtering properties in qualitative agreement with experiments. It should be easy to make 
this model more realistic by incorporating different radial functions of o and £ when there will be 
better constraints by measurements. 

The present results lead to several interesting perspectives for future work or extensions. First, as 
mentioned above, the simplified model could be enhanced by comparison with a more realistic model, 
for example based on three-dimensional reconstructions of extracellular space. The simplified model 
could be adjusted so that it fits as closely as possible the behavior of the more realistic model, yielding 
more optimal expressions of the radial profiles of conductivity and permittivity. A second possible 
direction is the "reverse" problem of estimating neuronal activity based on LFP measurements. By 
using data on the spatial and temporal variations of multisite LFPs and multi-unit activity, it should be 
possible to estimate what are the respective contributions of the natural frequency-filtering properties 
of extracellular space and the spatial coherence of neuronal sources in the different frequency com- 
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ponents. For example, it was shown that a consistent relation between LFP and cell firing extends to 
large cortical distances (>7 mm) for slow-waves but not for fast oscillations in the gamma (20-60 Hz) 
frequency range (Destexhe et al., 1999). A model of LFP is needed to evaluate whether this effect 
is really due to differences in the coherence of neuronal firing (as the single-unit data indicates), or 
if a large part could be explained by the low-pass filtering properties of the extracellular medium. A 
combination of experimental recordings and computational models will be needed to understand how 
neuronal activity translates into extracellular field potentials and vice-versa. 
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Appendix 1: Extracellular potential in non-homogeneous media with spherical 
symmetry 

This Appendix refers to Section Stationary currents in spherically-symmetric non-homogeneous medium: 
Recall a = o(r), £ =const. Then Gauss' law becomes 

V-E = -. (38) 
£ 

The law of charge conservation becomes 

= V-(oE) = gV-E + E-(Vo) 

= °p + E-(Va). (39) 
Using Ohm's law and assuming o^O, this becomes 



at £ o £ 



We make the assumption that the current density j is stationary, i.e. it does not explicitly depend 
on time. Also, by assumption a and £ are time independent. Hence denoting a = | and b = j • VlogG, 
Eq.|40|takes the form 

^- = -ap-b = 0. (41) 
dt 

Consequently, the induced charge density (at steady-state) is given by 

b 

P = — , (42) 
a 

p = -£E- VlogG. (43) 
It shows that the net charge density is different from zero. 

Electric field for spherical current source. 

At the radius r = R, the electric field is given by the current 

j\ r=R = o(R)E\ r=R . (44) 
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The total current passing through a sphere (surface S, radius r) is given by 



/= J dS-j = S(r)j(r) 



Thus the current density as function of radius r behaves as 

j(r) = j(r) e r = e r = e r 



S(r) 4nr 2 



Ohm's law implies for the electric field 



E(r)= ^ j(r)= 4^4w e ' 



Electric potential. 

Because the electric field is radially symmetric, so is also the potential, which obeys 

E(r) = -lv(r) . 
dr 



dr' 



Its solution is obtained from Eq.HTl 

V(r) = - I dr' E(r 



Anr' 2 o(r') ' 



(45) 



(46) 



(47) 



(48) 



(49) 



An independent check of this solution can be obtained by considering the law of charge conservation. 

(50) 

Because of the time-independent charge density, this becomes 



V.(o*> = -* 



V ■ (gW) = (Vo) • W + oAV = . 



Radial symmetry implies the following differential equation in the variable r, 



(51) 



(d r <5(rMd r V(r))+<5(r) 



-d r V + d 2 r V 



= 0. 



(52) 



A straight forward calculation shows that the potential, given by Eq.|4_9]is the solution of this differ- 
ential equation for V — ► when r — ► °°. 
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Appendix 2: Fourier component of the extracellular field potential 

This Appendix refers to Section Time-varying currents in non-homogeneous medium: 
Starting from Gauss' law 

V-(eE)=p, (53) 

the inhomogeneity of £ implies 

E-(Ve)+£V-E = p . (54) 

Then the electric potential obeys 

-(W)-(Ve)-eAV = p. (55) 

We recall that the potential V and the charge density p are time-dependent, while the permittivity £ is 
not. We define the Fourier transform of time-dependent function f(t) via 

poo 

U= dt e' at f(t) . (56) 



Now we perform a Fourier transform with respect to time of the potential V and the charge density p 
to obtain an equation for the Fourier components at frequency oo, 

(Vy C0 )-(V£)+£A\/ C0 = -p a) . (57) 

Similarly, starting from the differential form of the law of charge conservation, 

9p 

dt 



V-(gE) = -^, (58) 



the inhomogeneity of a implies 



Then the potential obeys 



E-(Vg)+gV-E = . (59) 

at 



(W)-(Vo)+oAV= ^ . (60) 



dt 

From the Fourier transform of p follows that a component at frequency CO satisfies 



^rpco = topco • (61) 
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The last two equations imply that the Fourier component at frequency co of the potential obeys 

( Wo,) ■ (Vo) + oAV a = -Uop m . (62) 

Combining Eqs . 1571 and l62l yields 

(Wo) • (V(o + icoe)) + (a + icoeJAVo, = . (63) 

Appendix 3: Method to calculate the extracellular field potential from point 
current sources 

1. Compute the Fourier component (0 of the impedance 

Z (r) - -J— Hdr' -I [ q (*) + fme (*)] (64) 

where (0 = 2^/. This expression incorporates the values of the conductivity c(r) and permit- 
tivity e(r) as a function of the distance r. It is also assumed that V(°°) = 0. 

This quantity is computed for each frequency component a) of the spectrum, and for each extra- 
cellular distance r considered. It can be precalculated and stored in a matrix (Z[f] [r]). 

2. For each current source, compute the (complex) Fourier transform of the total membrane cur- 
rent, which we call here I^. 

3. For each current source, compute the Fourier component (0 of the extracellular potential: 

Vco(r) = Z M (r) /co . (65) 

4. For each current source, compute the extracellular potential by applying the (complex) inverse 
Fourier transform to Eq.l65l 



5. Finally, combine the contributions from all current sources to yield the extracellular potential at 
a given position x in the extracellular space. 
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Figure 1 : Radial variations of conductivity and permittivity can induce frequency-filtering properties. 
A. Scheme of the current source in radial symmetry. The current source is assumed to be spherical 
(continuous line; radius R). The conductivity and permittivity vary in radial symmetry according to 
the distance r from the center of the source. B. Conductivity a vs. radial distance r. Two cases are 
shown: (1) o(r)/o(R) = 1 + and (2) o(r)/o(R) = 1 - yfr^fr, where r = 0.2025 R (i?=l 

here). C. Permittivity £ vs. radial distance r. The two curves shown are: (3) E(r)/o(R) = 0.01 and (4) 
e(r)/c{R) = 0.01 [1 — a/ tq / r] . D-F. Real part (D), imaginary part (E) and norm (F) of the impedance 
Z co (r = 5R) vs. frequency /. Combining the profiles (1) and (3) in B-C leads to a high-pass filter 
(dashed line), whereas (2+3) gives low-pass characteristics (solid line). The combination (2+4) is 
such that c(r)/e(r) = const., in which case there is no frequency dependence (dotted line). 

Figure 2: Frequency-filtering properties obtained by a localized drop in conductivity. 
A. Profile of conductivity vs. distance. The conductivity was described by o(r)/o(R) = 1 — 0.2 (r — 
6R)/R for 6R < r < UR, o(r)/o(R) = -1 +0.2 (r-6R)/R for 117? < r < 16R, and o(r)/a(R) = 1 
otherwise. B. Profile of permittivity. e(r)/o(R) was constant and equal to 0.01. C-E. Real part (C), 
imaginary part (D) and norm (E) of the impedance as a function of frequency /. Z co (r) is shown 
for different distances r away from the source. F. Attenuation of the impedance norm |Z ro (r)| with 
distance. The different curves correspond to three different frequencies. 
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Figure 3: Frequency-filtering properties obtained from a periodically varying conductivity. 
A. Oscillatory profile of conductivity vs. distance (o(r)/o(R) = 0.501 +0.5 *cos[2n(r — R)/2R\. B. 
Profile of permittivity (e(r)/c(i?) = 0.01). C-E. Real part (C), imaginary part (D) and norm (E) of the 
impedance Z m (r) vs. frequency /. The different curves are taken at different distances r outside of 
the current source. F. Attenuation of the impedance norm |Z(o(r)| with distance. The different curves 
indicate the attenuation obtained at different frequencies. 



Figure 4: Frequency-filtering properties obtained with exponential decrease of conductivity. 
A. Profile of conductivity. o(r)/o(R) decays exponentially according to o(r)/c(R) = Co + (1 — 
Co) exp[—(r — R) fX\, with a space constant X = 500 /urn. B. Profile of permittivity. e(r)/c(/?) was 
constant (0.01). C-E. Real part (C), imaginary part (D) and norm (E) of the impedance Z co (r) vs. fre- 
quency /. The different curves show the impedance calculated at different distances r. F. Attenuation 
of the impedance norm |Z ra (r) | with distance. The different curves indicate the attenuation obtained 
at different frequencies. 



Figure 5: Distance dependence of frequency-filtering properties 

A. Ratio of impedance at fast and slow frequencies (Qioo) represented as a function of distance r 
(units of R). The <2ioo ratios are represented for different profiles of conductivity. Drop: localized 
drop of conductivity (short dash; same parameters as in Fig . |2l» . Osc: oscillatory profile of conductivity 
(solid line; same parameters as in Fig.|3l the dotted line indicates a damped cosine oscillation). Exp: 
exponential decrease of conductivity (long dash; same parameters as in Fig.HJexcept R = 1> X=10R). 

B. Profiles of conductivity with exponential decay (same parameters as in Fig.0 space constants X 
indicated in pm). C. <2ioo ratios obtained for the conductivity profiles shown in B. 
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Figure 6: Frequency-filtered extracellular field potentials in a conductance-based model. 
A. Membrane potential of a single-compartment model containing voltage-dependent Na + and K + 
conductances and a glutamatergic synaptic conductance. The glutamatergic synapse was stimulated 
at t = 5 ms (arrow) and evoked an action potential. B. Total membrane current generated by this 
model. Negative currents correspond to Na + and glutamatergic conductances (inward currents), while 
positive currents correspond to K + conductances (outward currents). C. Power spectrum of the total 
current shown in B. D. Impedance at 500 /jm from the current source assuming a radial profile of con- 
ductivity and permittivity (as in Fig.©. E. Extracellular potential calculated at various distances from 
the source (5, 100, 500 and 1000 fjm). The frequency filtering properties can be seen by comparing 
the negative and positive deflections of the extracellular potential. The fast negative deflection almost 
disappeared at 1000 jum whereas the slow positive deflection was still present. 
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